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ABSTRACT 


As  air  traffic  density  increases,  effective  use  of  airspace  must 
include  consideration  of  weather.  Accurate  identification  of  turbulent 
volumes  is  of  paramount  importance  to  flight  safety.  The  advent  of 
Doppler  radar  has  made  it  possible  to  observe  wind  motion  in 
convective  clouds.  Over  a  number  of  years,  research  has  lead  to  the 
spectral  width  (standard  deviation)  of  the  Doppler  velocity 
measurements  as  an  indicator  of  turbulence.  In  this  paper  we  address 
the  hypothesis  that  turbulence  is  essentially  isotropic  in  convective 
systems,  and  therefore,  observations  of  turbulence  are  independent  of 
viewing  angle. 

Radar  observations  made  during  the  months  of  April,  May,  and  June 
in  1980-1985  were  scanned  to  locate  storms  amenable  to  analysis.  A 
number  of  cases  are  presented  in  which  a  dual-Doppler  network 
provided  the  essential  data.  Each  case  is  in  a  different  quadrant  with 
respect  to  the  Norman  Doppler.  Forty-four  horizontal  planes  were 
studied  from  six  different  storms.  The  results  of  four  of  these  storms 
are  presented.  Maximum  reflectivity  of  these  storms  range  from  51  to 
58  dBZ.  Altitudes  included  in  this  study  range  from  near  surface  to  7 
km.  At  these  altitudes  the  maximum  spectral  width  was  12  ms'1.  For 
these  cases,  involving  nearly  30,000  data  points,  70%  of  the  spectral 
width  observations  from  CIM  and  NRO  were  within  1  ms-1,  and  88%  had 
a  difference  of  2  ms'1  or  less.  These  results  indicate  that  the  use  of 
Doppler  radar  to  detect  turbulent  regions  has  a  high  probability  of 
success,  and  the  turbulent  regions  can  be  detected  independent  of  the 
direction  from  which  they  are  viewed. 
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1983  NSSL  Spring  Program  Observational  facilities  showing 
dual-Doppler  network 

Location  and  dates  of  analyzed  storms  relative  iu  the  mdar 
sites  at  Norman  OK  (NRO)  and  Page  Field,  Oklahoma  City,  Ok 
(CIM) 

NSSL  WSR-57  radar  scope  for  1628CST  17  May  1983. 
Elevation  angle  is  0.2  deg.  Range  marks  are  at  40km 
intervals.  Echo  contour  levels  are  at  lOdBZ  intervals 
starting  at  20dBZ.  Storm  studied  is  indicated  by  arrow. 

Section  of  the  computed  radar  reflectivity  (in  dBZ)  at  2km 
AGL.  17  May  1983  1628CST  as  seen  by  CIM.  Spacing  of  data 
points  is  1km. 

Correlation  coefficients  for  CIM-NRO  reflectivity  field  at 
2km  17  May  1983. 

Distribution  of  reflectivity  differences  (dBZ)  between  NRO 
and  CIM  for  all  altitudes  17  May  1983  1628  CST. 

Section  of  computed  CIM  spectral  widths  (ms-1)  at  2  km  17 
May  1983  1628  CST.  Data  are  at  1  km  intervals. 

Spectral  width  differences  (ms-1)  at  2  km  17  May  1983 
corresponding  to  area  in  Figures  4  and  7. 

Summary  of  spectral  width  differences  between  NRO  and 
CIM  for  all  altitudes  combined.  Shown  is  the  distribution  of 
the  8500  data  points  quantized  by  reflectivity  at  20-29dBZ, 
30-39dBZ  and  over  40dBZ.  In  addition  the  combined 
distribution  is  shown  with  a  heavy  black  line. 
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Figure  10  17  May  1983,  1628CST,  2km  level  with  hatched  areas 
showing  vertical  motion  centers  of  6ms-1  or  more,  and 
dotted  areas  corresponding  to  spectral  widths  regions  of 
6ms-1  or  more.  Solid  contour  lines  indicate  reflectivity 
levels  (dBZ). 

Figure  11  NSSL  WSR-57  weather  radar  scope  photograph  25  May  1983 
1707  CST.  Range  marks  at  40  km  intervals.  Reflectivity 
contours  at  lOdBZ  intervals  starting  at  20dBZ.  Arrow 
indicates  storm.  Numerous  small  echoes  in  photograph  are 
from  aircraft  transponders  showing  location  of  aircraft. 

Figure  12  Cross  correlation  coefficients  between  NRO  and  C1M 
reflectivity  fields  for  5km  25  May  1983  1707  CST. 

Figure  13  Section  of  computed  NRO  Doppler  radar  spectral  widths 

(ms-1)  for  5  km  25  May  1983  1707  CST.  Data  are  at  1  km 
intervals. 

Figure  14  Cross  correlation  coefficients  between  NRO  and  CIM 

spectral  width  fields  for  the  same  scan  as  in  Figure  12. 

Figure  15  Section  of  spectral  width  differences  (ms-1)  between  NRO 
and  CIM  at  5  km  25  May  1983  1707  CST. 

Figure  16  Summary  of  spectral  width  differences  (ms-1)  for  all 

altitudes  in  storm  25  May  1983  1707  CST.  Heavy  dark  line 
is  the  integrated  value  for  all  reflectiveness.  (11,250  data 
points). 

Figure  17  Similar  to  figure  11  except  for  27  may  1983  1740  CST. 

Figure  18  Cross  correlation  coefficients  between  NRO  and  CIM 
reflectivity  fields  for  2  km  27  May  1983. 

Figure  19  CIM  reflectivities  (dBZ)  at  2  km  27  May  1983  1740  CST. 
Data  are  at  1  km  spacing. 
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Figure  20  Computed  2  km  spectral  width  field  (ms-1)  observed  by  CIM 
corresponding  to  figure  19. 

Figure  21  Spectral  width  differences  (ms-1)  between  NRO  and  CIM  at 
2  km  27  May  1983  1740  CST. 

Figure  22  Summary  of  spectral  width  differences  for  all  altitudes. 

Similar  to  figure  16  except  for  27  May  1983.  (720  data 
points). 

Figure  23  WSR-57  radar  scope  photograph  similar  to  figure  11  except 
for  28  May  1983  2008  CST. 

Figure  24  Cross  correlation  coefficients  between  NRO  and  CIM 

reflectivity  fields  at  3  km  for  28  May  1983,  2008  CST. 

Figure  25  Section  of  the  CIM  spectral  width  field  at  3km  28  May  1983 
2008  CST.  Data  are  at  1  km  intervals. 

Figure  26  Spectral  width  differences  (ms-1)  between  NRO  and  CIM 
corresponding  to  data  in  figure  25. 

Figure  27  Summary  of  spectral  width  differences  (ms-1)  for  all 

altitudes  28  May  1983  2008  CST  occurring  at  specific  dBZ 
intervals.  Heavy  dark  line  is  the  distribution  for  all 
reflectivities.  (9,450  data  points). 

Figure  28  Composite  of  spectral  width  differences  for  the  four 
storms.  (29,850  data  points). 
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TURBULENCE  SPECTRAL  WIDTHS  VIEW  ANGLE  INDEPENDENCE 
AS  OBSERVED  BY  DOPPLER  RADAR 

J.  T.  Lee  and  K.  Thomas 

Cooperative  Institute  for  Mesoscaie  Meteorological  Studies 
University  of  Oklahoma 
401  East  Boyd 
Norman,  OK  73019 


1.  INTRODUCTION 

Aircrafts  usually  avoid  large  areas  of  airspace  occupied  by 
thunderstorms  resulting  in  service  disruption  and  increased  fuel 
consumption.  At  times  this  means  missed  connections  and/or  cancelled 
flights.  The  economic  considerations  may  influence  decisions  as  to  flight 
paths  which  may  turn  out  to  be  a  disadvantage  to  aircraft,  crew,  and 
passengers.  The  advent  of  Doppler  radar  offers  an  important  new 
capability  for  addressing  this  problem. 

Over  a  number  of  years,  cooperative  programs  between  the  Federal 
Aviation  Administration  (FAA),  the  United  Stated  Air  Force  (USAF)  through 
the  Air  Force  Cambridge  Research  Laboratory  (AFCRL)  and  the  Aeronautical 
Systems  Division  (ASD),  the  National  Aeronautics  and  Space 
Administration  (NASA),  universities,  including  Massachusetts  Institute  of 
Technology’s  Lincoln  Laboratory,  the  National  Center  for  Atmopsheric 
Research  (NCAR)  and  the  National  Oceanic  and  Atmospheric 
Administration’s  National  Severe  Storms  Laboratory  (NSSL)  have  been 
conducted  to  observe  thunderstorm  turbulence  with  aircraft  and  Doppler 
radar. 

The  results  of  these  programs  have  been  published  (Lee,  1967;  Burnham 
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and  Lee,  1969;  Lee  and  Carpenter,  1979;  Doviak  and  Lee,  1985).  It  was 
shown  that  turbulence  experienced  by  aircraft  could  be  related  to  the 
spectral  width  of  the  radial  velocity  measured  by  the  Doppler  radar 
(Donaldson  and  Wexier,  1969,  Lee,  1977;  Labitt,  1981;  Bohne,  1982;  Lee. 

Y.,  et  al,  1985;  and  others).  It  was  also  found  that  spectral  widths 
corresponding  to  turbulence  greater  than  moderate  were  observed  in  only 
about  30%  of  the  volume  of  a  tornadic  storm  -  that  is  70%  of  the  storm 
may  have  only  light-to-moderate  turbulence  (Doviak  et  al.,  1978).  If  the 
intensity  and  location  of  turbulence  can  be  accurately  and  routinely 
detected  by  ground  based  Doppler  radar,  aircraft  can  be  safely  routed 
around  the  hazardous  regions  (Zrnic'  and  Lee,  1982;  Mahapatra  and  Lee, 
1984). 

The  viewing  angle  independence  is  thus  important.  Spectral  width  is  a 
scaler  parameter  which,  in  isotropic  turbulence,  is  independent  of  the 
radar's  viewing  angle.  One  purpose  of  this  study  is  to  examine  to  what 
extent  the  turbulence  in  the  storms  is  nearly  isotropic.  This  is 
accomplished  by  comparing  spectral  width  measurements  made  by  two 
Doppler  radars  separated  by  tens  of  kilometers,  and  viewing  the  same 
storm  from  a  different  angle. 

2  SPECTRAL  WIDTH 

Atmospheric  turbulence  in  convective  clouds  is  a  fluctuation  of  the 
wind  around  a  mean  value.  These  fluctuations  are  caused  by  mechanical 
and/or  thermal  forces  acting  independently  or  in  concert.  In  a  turbulent 
flow,  turbulent  energy  is  assumed  to  cascade  from  the  longer  wavelengths 
(large  eddy  sizes)  to  the  shorter  wavelengths  (smaller  eddy  sizes)  in  a 
dissipation  process  (MacCready,  1964;  Frisch  and  Strauch,  1976;  Gage, 
1979;  Lilly  and  Peterson,  1983;  Brewster  and  Zrnic',  1986). 
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Doppler  radar  can  measure  these  variations  in  the  observable  wind 
field.  These  variations  within  a  resolution  volume  will  result  in  an 
increase  in  the  spectral  width  of  the  Doppler-estimated  mean  wind. 

The  Doppler  radar’s  spectral  width  estimation  is  a  function  of  radar 
system  parameters  (Zrnic’,  1977  and  1979)  and  meteorological  factors. 
There  are  four  major  mechanisms  which  can  effectively  broaden  the 
spectral  widths  --  that  is,  increase  its  numerical  value.  These  are: 
antenna  angular  velocity  rotation  (r),  differences  in  drop  sizes  (d),  and 
therefore,  differences  in  fall  speeds,  radial  velocity  shears  (s),  and 
turbulence  (t).  This  relationship  is  expressed  as: 

s2(total)  =  s2(r)  +  s2(d)  +  s2(s)  +  s2(t), 
where  s2(total)  is  the  total  spectral  width  squared.  The  beam  broadening 
terms,  as  well  as  the  effect  of  drop  size  distributions,  have  been 
evaluated  by  a  number  of  researchers  (Nathanson,  1969;  Lhermitte,  1963; 
Doviak  et  al.,  1979).  These  studies  indicate  that  the  contributions,  by 
beam  broadening  (antenna  rotation)  and  different  drop  size  distributions 
can  effectively  be  neglected.  However,  there  exist  certain  conditions  of 
strong  shear  that  will  cause  the  shear  term  to  be  significant. 

In  this  study,  we  will  assume  that  the  turbulence  is  homogeneous 
within  the  resolution  volume,  defined  as  that  volume  within  the  half¬ 
power  interval  of  the  radiated  beam  and  the  depth  of  the  individual 
sampling  interval.  For  further  discussion  of  turbulence  and  resolution 
volume  see  Doviak  and  Zrnic'  (1984)  and  istok  and  Doviak  (1936). 

To  utilize  this  information  and  to  test  if  the  turbulence  is  isotropic  or 
nearly  isotropic,  and  therefore  independent  of  viewing  angle,  a  series  of 
observations  were  made  using  a  configuration  of  two  Doppler  radars.  This 
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dual-Doppler  network  is  shown  in  Figure  1.  The  radars  located  at  NSSL  in 
Norman  (NRO)  and  at  Page  Field  (formerly  Cimarron  Airport)  (CIM)  are 
separated  by  41.473  km.  CIM  is  on  the  310.08  deg.  radial  (true)  from  NRO. 
The  individual  characteristics  of  these  radars  are  given  in  Table  1. 

The  Dopper  data  are  processed  as  detailed  in  Appendix  A.  Generally  the 
data  are  assigned  to  grid  points  spaced  at  1  km  intervals  (Brown  et  al, 
1981),  both  in  the  vertical  and  the  horizontal.  A  Cressman  spherical 
objective  analysis  scheme  is  used  with  a  1  or  1.5  km  radius  of  influence, 
depending  on  the  Doppler  radar  scanning  strategy  on  a  particular  day.  A 
minimum  of  two  data  points  from  the  same  radar  is  required  for  inclusion. 
Areas  of  insufficient  data  are  deleted.  Data  in  regions  where  the 
reflectivity  is  less  than  20  dBZ  are  not  used  in  the  final  analysis  because 
of  possible  noise  contamination. 

The  following  cases  are  presented  to  illustrate  the  independence  of 
spectral  width  from  viewing  angle  restrictions. 

3.  TEST  CASES 

A  number  of  cases  were  studied.  Four  are  presented  which  illustrate 
various  viewing  angles,  different  ranges,  and  several  elevation  angles. 
Each  case  is  in  a  different  quadrant  from  the  Norman  radar.  Figure  2 
shows  the  relative  position  of  these  storms  and  the  Doppler  radars. 

3.i  May.jz  iaaa 

On  May  17,  1983,  the  synoptic  situation  supports  a  strong 
southwesterly  flow  aloft  and  a  dryline  in  western  Oklahoma. 
Thunderstorms  (some  tornadic)  first  develop  in  southwestern  Kansas. 
Others  soon  develop  in  northwest  Oklahoma.  These  thunderstorms  move 
eastward  and  new  cells  develop  southward,  forming  a  NNE-SSW  line.  By 
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1628  CST,  this  line  moves  to  about  40  km  northwest  of  Norman. 


Figure  3  is  a  photograph  of  the  10-cm  WSR-57  weather  radar  PPI 
display.  This  radar  is  located  at  Norman  adjacent  to  the  Doppler  radar. 

The  storm  covered  by  the  dual-Doppler  observations  is  located  at  A,  at 
265°/50  km  from  NRO,  and  208°/38  km  from  CIM.  Thus  the  beams 
intersect  at  an  angle  of  41°.  Beam  width  at  this  range  is  0  7  km  for  NRO 
and  0.5  km  for  CIM.  Maximum  reflectivity  factor  of  the  storm  at  this  time 
is  54  dBZ.  The  maximum  total  spectral  width  (a)  is  10  ms'1. 

The  dual-Doppler  data  are  processed  and  grided.  Six  horizontal  planes  - 
surface  (lowest  elevation  available)  1,  2,  3,  4,  and  5  km  above  ground, 
with  approximately  10,000  data  points  -  are  examined.  Initially  the  NRO 
and  CIM  reflectivities  are  compared  and  cross-correlated  to  establish  if 
there  are  any  range  or  azimuth  discrepancies  besides  the  ones  accounted 
for  by  calibrations.  For  May  17th,  a  maximum  correlation  coefficient  of 
97%  is  observed  for  a  0,0  lag.  In  Figure  4  we  show  the  radar  reflectivity 
at  2  km,  as  seen  by  the  Cimarron  radar.  Figure  5  is  the  corresponding 
correlation  coefficient  as  the  data  are  lagged  in  both  the  x  and  y 
directions.  Results  are  similar  at  other  altitudes,  therefore  no  range  or 
azimuth  correction  is  applied  to  these  data.  The  distribution  of  the 
differences  between  the  corresponding  reflectivity  values  for  all  data 
points  is  given  in  Figure  6,  and  we  shall  thus  assume  that  the  radars  are 
observing  the  same  volumes. 

The  corresponding  analysis  is  then  performed  on  the  spectral  widths. 
Figure  7  is  an  example  of  the  spectral  width  field  for  the  2  km  level,  as 
observed  by  the  CIM  radar.  In  the  portion  of  the  storm  displayed  in  Figure 
7,  the  spectral  width  varies  from  2  to  9  ms'1,  with  4  ms'1  being  the 
dominant  value.  Note  that  in  the  region  of  50  dBZ  reflectivity  values  in 
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Figure  4  (coordinates  x  =»  -19.14  and  y  -27.7),  the  corresponding  spectral 
widths  are  3-4  ms*1,  indicating  light-to-moderate  turbulence  (Donaldson 
and  Wexler,  1969;  Lee,  1979).  Some  of  the  large  spectral  widths  of  9  ms'1 
at  x  =  -24.14  and  y  »  45.7  are  most  probably  due  to  artifacts  at  low  signal 

strength.  Figure  8  shows  the  corresponding  differences  between  the  CIM 

1 

and  NRO  observations.  In  general,  these  differences  are  within  ±2  ms*1.  A 
few  4  ms'1  values  at  x  =  -57  ,  y  =  +7  and  at  x  «  -57,  y  »  -16  are  at  the 
edges  of  the  storm  and  are  the  result  of  contamination  of  the  NRO  radar 
data  by  noisy  signal.  Similar  comparisons  are  made  for  the  other  five 
horizontal  planes.  These  results  are  combined  and  summarized  in  Figure  9. 
The  heavy  black  line  is  for  all  the  data.  Note  that  more  than  70%  of  the 
observations  are  within  1  ms*1,  and  86%  are  within  2  ms*1.  In  this  figure 

we  also  show  the  distribution  of  the  differences  classified  according  to 

I 

reflectivity.  From  this  figure,  there  appears  to  be  very  little  reflectivity 

I 

dependence.  Thus  is  this  storm,  with  reflectivities  ranging  up  to  54  dBZ 
and  at  altitudes  up  to  5  km,  the  observed  spectral  widths  obtained  by  NRO 
and  CIM  are  essentially  the  same  when  compared  point-by-point. 

The  vertical  motion  in  the  storm  is  calculated  from  the  dual-Doppler 
data  following  the  method  of  Brown,  et  al.  (1981).  Figure  10  shows  the 
relative  position  of  the  radar  reflectivity  contours,  areas  of  significant 
spectral  widths,  and  the  vertical  motion  fields.  Particular  notice  should 
be  made  that  the  areas  of  maximum  reflectivities  are  not  necessarily 
colocated  with  the  areas  of  maximum  spectral  widths.  Indications  are 
that  the  turbulence  is  more  closely  associated  with  the  edge  of  updrafts 
and  downdrafts.  Frisch  and  Strauch  (1976),  and  Knupp  and  Cotton  (1982) 
showed  similar  patterns.  The  future  of  forecasting  turbulent  areas  may 


well  depend  on  the  success  ir.  detailing  these  vertical  motions  within  the 

storm. 

3.2  May  25.  1.983 

On  25  May  1933,  moist,  unstable  air  is  brought  into  central  Oklahoma 
by  low-level  southeasterly  winds.  Surface  dew  points  are  near  14  °C.  A 
north  westerly  wind  at  500  mb  guides  a  cold  front  from  northern 
Oklahoma  through  Norman  and  Tulsa  by  mid-afternoon.  Thunderstorms 
form  on  the  front  east  of  Norman  and  the  main  center  of  activity  is  120 
km  southeast  of  Norman  at  1700  CST.  The  storm  covered  by  the  dual- 
Doppler  observations  is  at  138°/103  km  from  NRO,  and  135°/144  km  from 
CIM.  The  radar  beams  intersect  at  an  angle  of  3°.  The  NRO  beamwidth  at 
the  storm  is  1.4  km  and  corresponding  value  for  CIM  is  1.9  km.  In  this 
case,  CIM  is  looking  almost  directly  over  the  shoulder  of  NRO.  Figure  11  is 
a  photograph  of  the  WSR-57  radar  scope  at  0°  elevation  angle  at  1707  CST. 
Maximum  reflectivity  of  the  storm  is  58  dBZ  and  the  spectral  widths  vary 
from  2  to  10  ms-1.  Due  to  the  greater  distance  from  CIM  on  May  17th,  the 
lowest  usable  plane  is  at  3  km.  Data  for  3,  4,  5,  6,  and  7  km  are  used  in 
the  analysis.  Approximately  12,000  data  points  are  in  the  grid. 

Again  we  test  the  radars  for  range  and  azimuth  consistency  and  obtain 
a  maximum  correlation  coefficient  of  97%  at  0,0  lag  (Figure  12). 

Figure  13  is  a  portion  of  the  spectral  width  field  at  5  km,  as  observed 
by  CIM.  For  this  case  the  widths  are  generally  larger  (7,  8,  and  9  ms'1) 
than  in  the  previous  case.  This  is  to  be  expected  since  the  storm  as  a 
whole  is  more  intense  than  the  one  on  May  17th.  Figure  14  shows  the 
cross  correlation  coefficients  for  the  spectral  widths  for  CIM  and  NRO  at 
5  km.  The  69%  value  at  0,0  lag  is  not  as  large  as  for  the  reflectivity,  but 
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is  stilt  acceptable.  Figure  15  is  the  corresponding  difference  field  (NRO- 
CIM).  Here  again  we  find  good  correspondence,  with  82%  of  the  • 
differences  within  ±  1  ms*1  and  96%  within  2  ms*1  or  less. 

This  is  even  better  than  May  17th  even  though  the  CIM  beamwideth  is 
about  0,5  km  more  than  NRO's  due  to  the  greater  distance  to  the  storm 
from  CIM.  The  high  correspondence  is  very  encouraging. 


For  27  May  1983,  northwesterly  flow  at  500  mb  continues  to  dominate 
the  state.  At  the  surface  the  air  is  moist  (dew  point  about  *i5°C)  and 
unstable.  There  is  a  north-south  quasi-stationary  front  in  eastern 
Oklahoma  with  strong  southerly  winds  to  the  west  of  the  front,  including 
the  Oklahoma  City  area.  In  late  afternoon  several  thunderstorms  form 
northeast  of  Oklahoma  City.  Dual-Doppler  coverage  at  1740  CST  is  1 
centered  on  a  storm  at  40°  and  64  km  from  Norman.  From  CIM  the  azimuth 


and  range  is  70°  and  76  km.  The  NRO  beamwidth  near  the  storm  is  0.85  km 
and  CIM's  is  1.02  km.  These  beams  intersect  at  an  angle  of  32°.  Maximum 

reflectivity  of  the  storm  is  51  dBZ  and  3/4  inch  hail  is  reported  to  be 

;  j 

associated  with  this  storm.  Maximum  spectral  width  observed  is  8!  ms*1. 

| 

Figure  17  shows  the  WSR-57  radar  display  (0°  tilt  angle)  for  17^3  CST. 
Data  from  four  altitudes  --  near  tho  surface  (lowest  elevation  angle),  1, 

2,  and  3  km  --  are  available.  The  lowest  altitude  is  not  used  because  of 
contamination  of  the  CIM  data  by  returns  from  ground  objects  (ground 
clutter).  A  total  of  over  700  data  points  are  used.  This  storm  is  unique  in 
that  it  is  an  isolated  storm  and  the  possibility  of  second  trip  and  sidelobe 
contamination  is  minimized.  This  means  that  the  analysis  is 
straightforward  and  very  little  editing  is  required. 

Figure  18  illustrates  the  correlation  between  the  NRO  and  CIM 
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reflectivity  fields  at  an  altitude  of  2  km.  Once  again  the  highest 
correlation  is  at  0,0  lag,  indicating  that  there  has  been  no  detectable  shift 
in  range  or  azimuth  of  the  radar  data.  The  corresponding  2  km  CIM 
reflectivity  field  is  shown  in  Figure  19.  This  storm  is  basically  one  major 
cell  located  at  x  =  68,  y  =  28.  Figure  20  is  the  spectral  width  field  with  3 
to  4  ms-1  being  the  predominant  values,  except  near  the  northwest  edge  of 
the  storm  where  a  few  6  ms-1  are  found.  As  seen  in  Figure  21,  which 
shows  the  difference  in  the  spectral  widths  values  between  the  CIM  and 
NRO  observations,  this  northwest  area  also  has  the  largest  difference. 

Most  other  areas  have  values  near  0  ms'1.  The  graph  (Figure  22)  shows 
this  good  agreement  between  the  CIM  and  NRO  widths,  with  nearly  80%  of 
the  observational  differences  1  ms-1  or  less  and  98%  2  ms*1  or  less. 

3.4  Mav  28.  1983 

On  28  May  1983,  the  500  mb  flow  is  from  the  northwest.  Moist, 
unstable  air  is  moving  northward  over  the  state  during  the  day,  while  the 
remnants  of  a  mesoscale  convective  complex  (MCC;  large  area  of 
thunderstorms)  is  moving  into  western  Oklahoma  from  eastern  New  Mexico 
and  eastern  Colorado.  These  storms,  on  an  east  northeast/west  southwest 
line,  increased  in  intensity  and  by  1408  CST  the  leading  edge  was  about 
120  km  northwest  of  Norman.  The  storm  selected  for  analysis  is  317°  at 
118  km  from  NRO,  and  320°  at  78  km  from  CIM.  The  beams  intersect 
making  an  angle  of  3°.  In  this  case,  NRO  is  looking  over  the  shoulder  of 
CIM  --  the  reverse  of  the  May  25th  case.  The  NRO  radar  beamwidth  at  the 
storm  is  1.6  km  and  CIM's  is  1.0  km.  Maximum  reflectivity  of  the  storm  is 
55  dBZ.  Maximum  spectral  width  is  12  ms1.  Figure  23  is  the  2008  CST 
for  28  May  WSR-57  radar  scope  at  0°  elevation  angle.  Seven  altitude 
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section^  (1.  2,  3,  4,  5,  6,  and  7  km)  are  analyzed. 

in  F:n.  24  are  the  cross  correlation  coefficients  for  the  reflectivities 
and  in  Figs.  25  and  26  are  the  CiM  spectra!  widths  and  the  NRO-CIM 
spectial  width  differences  at  3  km.  Once  again  we  note  in  Figuie  24  that 
the  maximum  correlation  coefficient  is  near  the  0,0  lag,  but  a  slight  drift 
in  the  radar  range/azimuth  stability  begins  to  show  itself.  This 
particular  system  is  contaminated  by  large  areas  of  second  trip  data,  but 
the  storm  is  of  such  dimension  that  over  9,000  data  points  are  available 
for  comparison.  The  section  of  the  storm's  spectral  width  observations, 
shown  in  Figure  25,  has  values  of  1  to  9  ms'1  --  probably  the  largest  range 
of  values  for  any  of  the  cases.  We  note  that  the  spectral  width 
differences  in  Figure  26  indicate  no  particular  bias  for  either  radar,  but 
the  differences  are  greater  than  in  the  previous  cases,  with  only  58% 
within  1  ms'1  or  less  and  78%  within  2  ms'1  (Figure  27). 

This  larger  difference  can  be  ascribed  to  the  greater  range  to  the  storm 
(and  thus  the  greater  beamwidth)  and  to  the  fact  that  there  is  so  much 
range-overlaid  (second  trip  data)  which  was  not  flagged  by  the  computer 
program  (Brown,  et  al.,  1981).  This  increase  in  spectral  width  differences 
at  long  range  here  and  on  May  25th  suggests  that  a  possible  range 
limitation  may  have  to  be  determined. 

4.  SUMMARY 

Whereas  only  a  small  portion  of  the  analyzed  data  fields  have  been 
presented,  the  summary  graphs  (Figures  9,  16,  22  and  27)  show  that  there 
is  strong  evidence  that  turbulence  (spectral  widths),  as  detected  by 
ground-based  calibrated  radar,  is  independent  of  viewing  angle. 

Figure  28  is  a  composite  of  all  the  data  for  these  cases.  It 
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encompasses  altitudes  from  near  ground  level  to  7  km  and  ranges  to 
storms  from  35  km  to  120  km.  These  nearly  30,000  comparisons  show 
that  70%  of  the  NRO  and  CIM  spectral  width  observations  have  a  difference 
of  1  ms*1  or  less  and  88%  have  a  difference  of  2  ms'1  or  less.  The  mean  of 
the  distribution  is  -0.05  ms'1,  and  the  standard  deviation  is  2.66  ms"1. 

At  ranges  of  these  storms  with  reflectivities  more  than  20  dBZ,  the 
signal-to-receiver  noise  ratios  are  larger  than  20  dB.  This,  together  with 

a  number  of  samples  processed  to  obtain  ov’s,  implies  that  errors  in 
individual  spectral  width  estimates  from  either  radar  are  less  than 
0.9ms'1.  Because  at  least  two  estimates  from  each  radar  are  averaged  in 
the  process  of  interpolating  to  a  f  x  1  km  grid,  the  statistical  errors 
should  be  reduced  to  less  than  0.5  ms*1.  Thus  the  small  1  to  2  ms'1 
differences  are  real  and  are  most  likely  due  to  the  short  observation 
period  (0.1  s)  and  slight  inconsistencies  in  the  timing  of  the  two 
observations.  There  simply  may  not  be  sufficient  time  for  a  single 
recording  of  an  isotropic  turbulent  eddy  to  perfectly  resolve  this  property. 
Thus  based  on  the  data,  we  conclude  that  there  prevails  —  at  least  in  two 
dimensions  and  in  moderate-to-severe  storms  --  isotropy  of  turbulent 
eddies  with  scales  less  than  1  km.  This  is  important  because  it  means 
that  irrespective  of  the  location  of  the  turbulence,  if  there  are  scatters 
the  Doppler  radars  should  be  able  to  delineate  the  areas.  This  also 
suggests  that  enhanced  ability  to  detect  gust  fronts  and  downbursts 
outflows,  which  may  be  turbulent,  should  be  possible  by  the  inclusion  of 
the  spectral  width  in  the  algorithms,  since  outflows  in  themselves  may  be 
tangent  to  the  Doppler  beam  and  only  detected  by  the  spectral  width. 

The  next  step  after  detection  is  the  problem  of  forecasting  turbulent 
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APPENDIX 


DATA  ANALYSIS 

The  data  used  for  this  project  were  initially  processed  on  the  National 
Bureau  of  Standards'  CYBER  855.  Upon  completion  of  the  objective 
analyses,  the  data  were  sent  to  the  National  Severe  Storms  Laboratory's 
VAX  11/780,  where  the  remainder  of  the  analyses  were  completed. 

Doppler  data  were  processed  using  the  techniques  described  by  Brown  et 
al.  (1981). 

Prior  to  the  objective  analysis,  data  for  both  the  Norman  Doppler  (NRO) 
and  the  Cimarron  Doppler  (CIM)  were  edited  using  a  local  environmental 
check  algorithm  (Eilts  and  Doviak,  1986)  and  by  visual  inspection.  In  the 
local  environmental  check,  data  points  are  compared  to  the  surrounding 
data  with  questionable  data  points  automatically  removed. 

The  data  were  analyzed  using  the  Cressman  (1959)  weighting  scheme. 
The  radius  of  influence  used  was  dependent  of  the  data  resolution.  The 
grid  spacing  is  1  km,  in  the  horizontal  and  vertical.  Even  with  a  fairly 
large  vertical  influence  radius  there  were  a  few  problems.  At  extreme 
distances  (140  to  180  km)  from  the  radar  site,  data  points  in  the  vertical 
may  be  separated  by  distances  greater  than  the  objective  analysis  grid 
spacing  of  1  km.  This  problem  occurred  with  the  May  25,  1983  data  set, 
especially  with  CIM.  the  May  27th  NRO  data  also  fit  into  this  category.  At 
these  large  distances,  the  interpolated  values  may  not  be  representative. 
Consequently,  it  was  necessary  to  delete  some  data  values  at  large 
distances. 

Before  any  comparisons  are  performed,  several  consistency  checks  are 
made.  Cross  correlations  between  both  reflectivity  and  spectral  width 
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using  the  analyzed  data  fields  from  the  two  radars  are  computed,  fhe 
reason  for  iesting  the  correlation  is  to  make  sure  that  general  features 
are  not  displaced  in  range  or  azimuth.  In  this  procedure,  the  data  grid  of 
one  of  the  radars  is  shifted  around  by  1  km  at  a  time,  in  both  the  x  and  y 
directions,  and  new  correlations  calculated  using  the  following  equation: 


Correlation  [0+dx,  0+dy]  (P)  =  - fcU&xY)  - 

[N(ZX2)  -  (EX)2]1'2  [N(Lg2)  -  (Ey)2]1'2 

where  dx  and  dy  are  displacements  in  the  x  and  y  directions  of  the  CIM  data 
matrix  relative  to  the  NRO  matrix.  Ideally,  the  best  correlations  should 
occur  without  any  shifting  of  the  data  (0,0  lag).  For  this  study,  data  from 
the  two  radars  correlated  quite  well  except  for  the  May  28,  1983  case 
where  a  slight  correction  had  to  be  included. 

The  next  step  involves  checking  the  data  for  consistency  (Brown,  et  al., 
1981).  In  most  of  the  1983  cases,  the  NRO  and  CIM  reflectivity  values  are 
in  agreement.  I  i 

I 

Histograms  are  plotted  to  determine  if  a  correction  is  needed.  Irj  all 

•  t 

cases,  the  NRO  and  CIM  indicates  similar  distribution  of  the  data,  ekcept 
that  the  data,  at  times,  contained  an  error  value  offset,  this  error  value 
is  calculated  and  the  value  added  to  the  NRO  data.  Although  the  correction 
is  added  to  the  NRO  values,  this  does  not  imply  that  the  NRO  data  are 
necessarily  in  error.  The  corrections  used  are  indicated  in  the  Table  2. 

After  these  corrections  are  done,  a  manual  check  is  performed  that 
detects  the  obvious  bad  data  points.  Difference  fields  of  spectral  width 
are  then  generated  and  large  differences  are  noted.  At  these  suspect 
points,  both  NRO  and  CIM  data  are  examined  and  any  bad  points  are  deleted. 
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After  this  procedure,  the 

data 

are  considered  ready  f 

and  the  analysis  cornpio; 

;  'i '•  i  t  ..  .J 

•:.s  app:.ed  io  i\RO  Dcppr 

Dav.*  n-’c-.. 

,  .  .  r  . 

n  fiction  Spectral  Width 

May  17 

+  2  4 

+0.5 

May  25 

+7  6 

+  1.7 

May  27 

-r7.5 

+0.0 

May  28 

i-9  7 

+  1  4 
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Observational  facilities  showing  dual-Doppler  network 
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Figure  2  Location  and  dates  of  analyzed  storms  relative  to  the  radar 
sites  at  Norman  OK  (NRO)  and  Page  Field,  Oklahoma  City,  OK 
(CIM) 
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Figure  3  NSSL  WSR-57  radar  scope  photgraph  for  1628CST  17  May 
1983.  Elevation  angle  is  0:2  deg.  Range  marks  are  at  40km 
intervals.  Echo  contour  levels  are  at  lOdBZ  intervals 
starting  at  20dBZ.  Storm  studied  is  indicated  by  arrow. 
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Figure  4  Section  of  the  computed  radar  reflectivity  (in  dBZ)  at  2km 
AGL  17  May  1983  1628CST  as  seen  by  CIM.  Spacing  of  data 
points  is  1km. 
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Figure  5  Correlation  coefficients  for  CIM-NRO  reflectivity  field  at 
2km  17  May  1983. 


24 


Figure  6  Distribution  of  reflectivity  differences  (dBZ)  between  NRO 
and  CIM  for  all  altitudes  17  May  1983  1628  CST. 
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Figure  7 


Section  of  computed  CIM  spectral  widths  (ms-1)  at  2  km  17 
May  1983  1628  CST.  Data  are  at  1  km  intervals. 
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Figure  8  Spectral  width  differences  (ms-1)  at  2  km  17  May  1983 
corresponding  to  area  in  Figures  4  and  7. 
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Figure  9  Summary  or  spectral  width  differences  between  NRO  and 

CiM  fcr  al:  altitudes  combined.  Shown  is  the  distribution  of 
the  3500  data  points  quantized  by  reflectivity  at  2Q-29dl32, 
30-39dB7  and  ove;  40dBZ  in  addition  the  combined 
aistributior  o  shown  with  a  heavy  black  line 


17  MAY  1983 


2  KM 


Figure  10  17  May  1983,  1628CST,  2km  level  with  hatched  areas 
showing  vertical  motion  centers  of  6ms-1  or  more,  and 
dotted  areas  corresponding  to  spectral  widths  regions  of 
6ms-1  or  more.  Solid  contour  lines  indicate  reflectivity 
levels  (dBZ). 


Figure  11  NSSL  WSR-57  weather  radar  scope  photograph  25  May  1983 
1707  CST.  Range  marks  at  40  km  intervals.  Reflectivity 
contours  at  lOdBZ  inle-vals  starting  at  20dBZ.  Arrow 
indicates  storm  Numerous  small  echoes  in  photograph  are 
from  aircraft  transponder.':  showing  locator  of  aircraft 


Figure  13  Section  of  computed  NRO  Doppler  radar  spectral  widths 
(ms-1)  for  5  km  25  May  1983  1707  CST.  Data  are  at  1  km 
intervals. 
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Figure  14 


Cross  correlation  coefficients  between  NRO 
spectral  width  fields  for  tho  same  scan  as 


and  CIM 
in  Figure 
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Figure  15  Section  of  spectral  width  differences  (ms-1)  between  NRO 
and  CIM  at  5  km  25  May  1933  1707  CST.  .  ' 
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PERCENTAGE  OF  DIFFERENCE 


Figure  16  Summary  of  spectral  width  differences  (ms-1)  for  ail 

altitudes  in  storm  25  May  1983  1707  CST.  Heavy  dark  line 
is  the  integrated  value  for  all  reflectiveness.  (11,250  data 


Figure  17  Similar  to  figure  11  except  for  27  may  1983  1740  CST. 


36 


Figure  18  Cross  correlation  coefficients  between  NRO  and  C!M 
reflectivity  fields  for  2  km  27  May  1983. 


r 


Mil  i  p»  i  rt«  i,uii  mom 

■  Miif’unr  ippm 


mmmi  hr  nomicpmi  n*pf 


*  )  j  » *  •■•Pvt  6*0uhi> 


/.  »•*  'ill 


51.41 

41  42  4)  4* 


s  i 

1U  51  V?  H  H  Si  V 


»iii>in  <««j  »»pp  <i4»<(toii 

4'.  .«*  fu.ftl  /*-.45 


**ts  cmwiD  t *c»  ? io 


•  l/iiBD  lilt 

i  tO  »Q.O  6l|t 


•,..*  ft  1  62  S'*  ti.  ftS  tft  6  T  »4  69  73  / 1  72  M  7*  75  74  77  77  7t  10 


JB.4  64 

57.6  65 

36.6  *4 
11. 4  65 


54.6  42 

11.6  41  17  17 

12.4  40  «1  21 


11.  4  5*  24  70  ?J  71  7 1  19  14 

5Q.4  SB  27  16  51  14  5  5  2»  74  24  1  6 


24.4  57 
21. 6  56 

27.4  5S 
24.1  14 

21.4  55 
24.  4  52 
21.6  51 

72.4  10 

21.4  41 

20.4  41 

14.4  47 

11.4  44 

17.4  45 

14.4  4* 
11.1  45 

14.4  42 

15.4  41 

17.4  40 

11.4  1* 

10.4  14 

4.4  If 

1.4  14 

7.4  15 

6.4  54 


76  47  jv  40  4 «.  67  14  W  2 4  14 


26  65  44  46  66  45  64  56  11  25 
14  U  47  46  46  45  44  44  )4  1)  21 
21  40  51  64  42  40  40  44  4!  It  24  If 

74  41  45  41  16  11  18  41  64  61  55  24  20  24  II 

26  40  4 1  5«  51  51  5*  57  4t  42  5*  57  12  25  22 

70  52  <1  It  54  57  )|  56  56  57  54  40  57  55  56  2*  21  72  22 

21  54  41  40  57  54  J4  56  54  57  56  57  1«  54  17  11  51  26  21 

II  JO  51  54  54  11  JO  24  51  36  55  55  16  57  57  15  55  12  21  It 

24  54  54.  5«  it  24  21  26  55  54  54  55  55  55  55  ’4  52  2*  24  II 

25  21  2  4  74  24  II  #  27  24  52  55  14  J5  54  55  5  1  10  2*  25  21 

21  2  1  .  5  .*4  14  24  2 •  50  5  1  5  1  )5  5  }  5  2  51  50  .*6  2  1  II 

21  71  21  24  2»  >1  51  55  52  51  50  50  27  24  20 

54  24  27  24  11  52  52  51  2*  21  25  72  11 

16  22  27  21  79  11  51  24  24  25  22  14  17 

II  22  21  27  21  JO  10  27  25  25  20  17 

11  22  24  25  25  26  25  25  21  22 

1  II  20  22 


5W45 


ll'.ll 

I  IfSIMCt  Kll 


70*  4 1  75*  45 

non  (1MMI9  I0MUI  MMI 


1674 

•MINIUM  MM 
(MINI  TIM 
illl  CINTM 
L61HVII 
IQNtltllM 

If  IM41N  MON  HIM 
MIMNCt  MOP  IMM 
1  tOOMINAM 
1  < QQR4  tN  41 ( 


1711  II  Ml 
1740  If  Ml 

15  51  17. 1  N 
II  II  11.4  M 

71.1  Ml 
54.4  IN 
14.61  IN 
IP. 60  IN 


1101*  N0M0N 
ROIIM  MON 
MVtNI  61 

III!  POINt  IMC1NI 
I  It  MCI  I  OP 
I  »mcrt0N 
I  6IMCM0N 
Nil  IP »  Of  4«l»  MU 
(041  Mil 

con  riNi 


0.0  »U 
0.0  Mitt 

1.0  IN 
1.0  «N 
1.0  IN 
0.400  ««  *51. 
•5711720. 
15.55.00. 


40*41* 


Figure  19  CIM  reflectivities  (dBZ)  at  2  km  27  May  1983  1740  CST. 
Data  are  at  1  km  spacing. 
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Figure  20  Computed  2  km  spectral  width  field  (ms-1)  observed  by  CIM 
corresponding  to  figure  19. 
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Figure  21  Spectral  width  differences  (ms-1)  between  NRO  and  CIM  at 
2  km  27  May  1983  1740  CST. 
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Figure  22  Summary  of  spectral  width  differences  for  all  altitudes. 

Similar  to  figure  16  except  for  27  May  1983.  (720  data 

points). 
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Figure  23  WSR-57  radar  scope  photograph  similar  to  figure  11  except 
for  28  May  1983  2008  CST. 


42 


<  COflMCIENi;  ATI  0!iT»rHtO 
*  ’•  t  H  HOP1 1  f  f  *  111  A 

i »;*•  mi  t ) h  i  r h m or 

n-OH  (  OMMAB I  MG  DATA  .1 

1  A r  T 0 D  c»  103 

■jC  ai  ro  nm>)  i s  '  .u 

•i:  !  H  f 

'  Of 

t  .  i)  k; 

►  K  .  ..Ml) 

i.COtJ  KM.ftf  'tp|(  T  2  V r  LY 

*1  .  <«.  to 

•.  r . 

SI. 

«;  j  n 

* 

*  * 

5 0  . 

34. 

w. 

49. 

M  . 

5  4. 

5  5. 

5  7. 

f»  7  , 

5«, 

r  *  *  fc  . 

*  * 

*>* 

*  .  s  .  . 

.  » 

•  x  < 

, .  . 

3  4 . 

59. 

44. 

51 . 

56. 

iO. 

6  t  . 

64. 

65. 

65. 

i:  *  * 4 

.  .. 

..." 

•«*» 

.... 

/ 

5 

7,9 . 

•  • 

52. 

53. 

a5. 

nS. 

72. 

73. 

7*. 

71  , 

■  \  "■ 

.... 

* 

*  * 

... 

’  46. 

•  *  * 

53. 

»«*» 

61. 

65. 

75. 

eo. 

c2. 

50. 

77. 

73. 

3  3. 

61  . 

69, 

76. 

33. 

37. 

86. 

•  • 

78. 

72. 

| 

'  6  6. 

88. 

. 

’ '  *  ‘  ! 

56. 

65. 

. .  •  • 

74. 

51  . 

37, 

85. 

<50. 

74. 

•  * 

69. 

7*43. 

.  . . 

. . .  • 

»  •  *  *  * 

■ 

* 

.... 

******  1' 

56. 

65. 

74. 

81  . 

35. 

84. 

81. 

76. 

69. 

.5. 

'  f  5t. 

i 

"  *  * 

'  “ 

it. 

63. 

70. 

77.' 

79. 

78. 

74. 

70. 

64. 

58. 

5  2. 

■  * 

*  ■ 

4  _  * 

•  .  *  • 

*  *.* 

•  •  • 

* 

5  4  . 

39. 

63. 

71  , 

)S. 

73. 

69. 

6  3. 

53. 

u. 

47. 

*  *  *  * 

■  * 

... 

*  * 

*  *,* . 

52. 

55. 

61. 

66. 

68. 

67. 

63. 

58. 

51. 

46. 

f  41* 

Figure  24  Cross  correlation  coefficients  between  NRO  and  CIM 

reflectivity  fields  at  3  km  for  28  May  1983,  2008  CST 
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Figure  25  Section  of  the  CIM  spectral  width  field  at  3km  28  May  1983 
2008  CST.  Data  are  at  1  km  intervals. 
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Figure  26  Spectral  width  differences  (ms-1)  between  NRO  and  CtM 
corresponding  to  data  in  figure  25. 
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Figure  27  Summary  of  spectral  width  differences  (ms-1)  for  all 

altitudes  28  May  1983  2008  CST  occurring  at  specific  dBZ 
intervals.  Heavy  dark  line  is  the  distribution  for  all 
reflectivities.  (9,450  data  points). 
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Figure  28  Composite  of  spectral  width  differences  for  the  four 
storms.  (29,850  data  points). 


47 


